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Abstract 

We continue to investigate the physical interpretation of the Riemann zeta 
function as a FZZT brane partition function associated with a matrix/gravity 
correspondence begun in arxiv:0708.0645. We derive the master matrix of the 
(2, 1) minimal and (3, 1) minimal matrix model. We use it's characteristic poly- 
nomial to understand why the zeros of the FZZT partition function, which is 
the Airy function, lie on the real axis. We also introduce an iterative proce- 
dure that can describe the Riemann S function as a deformed minimal model 
whose deformation parameters are related to a Konsevich integrand. Finally 
we discuss the relation of our work to other approaches to the Riemann S 
function including expansion in terms of Meixner-Pollaczck polynomials and 
Riemann-Hilbert problems. 

1 Introduction 

Recently it has been remarked that the gauge/gravity correspondence is more fas- 
cinating from a physics point of view than the Riemann hypothesis [1]. Ironically 
in our approach, perhaps they are manifestations of the same phenomena. 

In a previous paper we interpreted the Riemann H function as the FZZT partition 
function of a matrix model [2]. On the gauge or matrix side of the correspondence 
the FZZT brane partition function is the expectation value of the characteristic 
polynomial of a matrix M weighted by e~ v ^ M ' ) with V(M) the potential of a matrix 
model [3] [4] [5] [6] [7] . On the gravity side it is the exponentiated macroscopic loop 
[8]. The matrix model in question is of a complicated type, but seems related to the 
(p, 1) minimal matrix models solved in [9] whose FZZT partition functions were the 
Airy function and generalized Airy functions [10]. These (p, 1) matrix models are 
special cases of the two matrix model [11] [12]. Historical papers on matrix theories 
and the large N limit are [13] [14] [15]. 
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In this paper we construct master matrices for these theories and discuss the 
application to the zeros of the FZZT partition function. A master matrix is a 
particular large N matrix whose characteristic polynomial yields the FZZT partition 
function in the large N limit [16] [17]. This method allows us to understand in 
particular why the zeros of the Airy function are on the real axis. We also introduce 
an iterative method to define the E function that makes the relation to the matrix 
model clearer. 

This paper is organized as follows. In section 2 we introduce the concept of the 
master matrix. We derive the master matrix of the (2, 1) minimal matrix model 
whose characteristic polynomial yields the FZZT partition function which is the 
Airy function. We also derive the master matrix of the (3, 1) minimal matrix model 
whose characteristic polynomial is the generalized Airy function. Interpreting the 
E function as a FZZT brane partition function, we discuss an iterative procedure 
which to leading order yields the deformed (3, 1) matrix model, to next order the 
deformed (5, 1) matrix model and so on. The deformation parameters are related to 
the expansion of the integrand of an integral representation of the E function which 
is similar to the integral representation of the generalized Airy functions. In that 
case the integrand is a Konsevich integrand which is cubic for the (2, 1) minimal 
model, quartic for the (3, 1) minimal model, sixth order for the (5, 1) minimal model 
etc. In section 3 we discuss the relation of our approach to other approaches to the 
Riemann E function including the expansion of the E function in terms Meixner- 
Pollaczek polynomials and the definition of the Riemann E function as the solution 
to a Riemann-Hilbert problem. In section 4 we state the main conclusions of the 
paper. 

2 Master matrix 

2.1 Master matrix of the (2, 1) minimal model 

The (2, 1) minimal model is defined by the partition function: 



and g is the coupling constant. In this paper we define a master matrix associated 
with the model as a matrix whose characteristic polynomial is equal to the matrix 
integral: 




with: 



V{M) = -Tr(M 2 ) 




which is the FZZT partition function. 
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Of particular interest is what happens as one takes the large N limit, as the 
zeros of the FZZT partition function are related to the eigenvalues of the master 
matrix. The master matrix for the (2, 1) minimal model is given by: 
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Which for N = 8 is given by: 
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The FZZT partition function for the (2,1) minimal model was computed in [10] and 
is: 

N/2 



(f) H N (z/J-g) 



This coincides with the characteristic polynomial of the master matrix. For the case 
N = 8 this is: 

105 4 105 , 9 105 o 4 fi o 

-^-5 4 - — g 3 z 2 + —g 2 z 4 - Ugz 6 + z 8 

The master matrix (2.1) agrees with the master matrix of the Gaussian matrix model 
computed in [16] which is has the same partition function as the (2, 1) minimal model 
after integration over P. 

Because the master matrix is manifestly Hermitian it's eigenvalues are real. The 
large N limit of FZZT partion function corresponds to [9]: 

9^ 7? 

and leads to the Airy function Ai(z). This function is given by the contour integral: 

v ' J 2m 

Co 
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with contour Cq starting at infinity with argument — 7r/3 and ending at infinity with 
argument 7r/3. It has the series expansion: 

Mz) = E 3 2 / 37r r((re y )/3) s[ ^ n + ik/sKs 1 ' 8 *)" 

n=0 

The Airy function obeys the differential equation: 

Ai"{z) - zAi(z) = 

We plot the Airy function on the real line in Figure 1 and in the complex plane in 
Figure 2. 
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Figure 1: Plot of the of the Airy function on the real axis. 



-iog(|Ai(.T+ ;>■)!) 




Figure 2: Plot of minus the logarithm of the magnitude of the Airy function in the 
complex plane. The zeros are all located on the negative real axis. 

The Airy function has all it's zeros on the real axis and this is a manifestation 
of the Hermitian nature of the master matrix in (2.1). 
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2.2 Master matrix of (3, 1) minimal model 

The (3, 1) minimal model is defined by the partition function with matrix potential: 

V(M) = i(^Tr(M 2 + ^Tr(M 3 )) 

The master matrix of the (3, 1) minimal model is the matrix M with nonzero com- 
ponents: 



(i - l)(i - 2)S itj+2 + 3(i - l)S i>j+1 + gS i+1 j 



which is of the form: 
/ 
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For N = 8 this is given by: 
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The characteristic polynomial of this master matrix for g = 1/N is given by: 

21z e 



8085 945z 175z 2 105z 3 945z 4 7z 5 
4096 ~ ^56 ~ + 16 + ~~ 32 F 



+ 2 C 



and this corresponds to the FZZT partition function of the (3, 1) minimal model 
computed in [9]. 



Q N {z)=g N {d x ) N e- 9 ^ 



ia x s + 3 x 2_ xz) 



\x=0 



for N = 8. The expression for Qn(z) can be written using the residue theorem as: 



Qn(z) = (-9) N N\^- j 



dip 



N+l 



-V{tp)-tpz 
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After taking the large N limit [9] : 



9^ W 



1 + J?tjZ Z 



one obtains a generalized Airy function 3>(z) denned by the integral: 



w J 2m 

Co 

Similar integrals in the context of black holes were considered in [18]. The general- 
ized Airy function obeys the differential equation: 

$"'(z) + z$(s) = 

with solutions: 

*(z) = A( F 2 ({},{^},-^)) 

+ B(z 2 ( 0J P 2 ({},{^},-^))) 

+ ^( F 2 ({},{^},-^))) 

for constants A, B, and C where pFq is a generalized hypergeometric function. The 
contour Co was chosen in [9] so that one obtains a solution which is real for real z 
and decays without oscillation for large positive z. 

Modifying the contour to be along the imaginary axis we can define a modified 
generalized Airy function ^f(z) by: 



e — A 4> +i<pz d( p 

-oo 

with a series expansion given by: 



Mz) = ±yt?lr( 1 - + k - 



z 2k 



This modified generalized Airy function obeys the differential equation: 

- zV(x) = 

with solution: 

m = ^n^n, { \, \ } , f 4 ) - . 2 r(^) F 2 ( {} , { \, 3 - } , 
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We plot the magnitude of this modified generalized Airy function ^(z) on the real 
axis in Figure 3 and in the complex plane in Figure 4. This function is even and 
as we shall see in the next section it has similar characteristics to the Riemann H 
function 




Figure 3: Plot of the of the logarithm of the magnitude of the function on the 
real axis. The position of the zeros are given by the location of the spikes pointing 
down. 




Figure 4: Plot of minus the logarithm of the magnitude of the Vl/ function in the 
complex plane. The zeros are symmetrically located on the positive and negative 
real axis. The position of the zeros are given by the location of spikes pointing up. 



2.3 Riemann H function 

The Riemann 3 function is defined by: 

E(z) = C(iz + -)T(i- + -W-^n-^t-- - I) 



7 



It is even and can be expressed as an integral along the imaginary axis as: 

/oo 
-oo 

where: 

OO Q 

U(<j)) = - log(X; (vr 2 fc 4 e^ - t^^y-^) 
k=i 2 

This function plays the same role for the 3 function as the Konsevich potential <p 3 /3 
plays for the Airy function and </> 4 /4 for the <3? function [19] [20] [21]. We plot the 
function U{4>) m Figure 5. It is even and this leads to the fact that the H function 
is even. For small <j) one can develop an expansion: 



1/(0) 




Figure 5: Plot of the function U{4>). The relation U(<p) = U(—(p) leads to the fact 
that the S function is even. 



U(<f>) = 9.36345^ + 5.95896</> 4 - 2.15104/ + 0(<p 8 ) (2.2) 

which is probably why the (3, 1) minimal model modified FZZT partition function 
shares some of the characteristics of the H function. 
The 3 function itself can be expanded as [22] : 



n=0 

where 



2n 

(2n)! * 



a 2n = 4jd£(£- 1 / i f(£)(ho g e) 2n ) 



and 



i 



oo 



f(£) = £ foV* ~ ^)£ 1/2 e-^ 
q =i 1 



Thus like the \P function one can think of the H(z) function as an infinite order 
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polynomial expanded in even powers of z. We plot the magnitude of the S function 
on the real axis in Figure 6 and in the complex plane in Figure 7. The Riemann 
hypothesis is equivalent to the statement that the zeros of the E function lie on the 
real axis. 




Figure 6: Plot of the logarithm of the magnitude of the H function on the real axis. 
The zeros are symmetrically located on the positive and negative real axis. The 
position of the zeros are given by the location of spikes pointing down. 




Figure 7: Plot of minus the logarithm of the magnitude of the 3 function in the 
complex plane. The zeros are symmetrically located on the positive and negative 
real axis. The position of the zeros are given by the location of spikes pointing up. 

The Riemann S function does not obey a finite order differential equation. Nev- 
ertheless keeping the first two terms in the expansion for U(4>) one can derive the 
following approximate equation for small z: 

4(5.95896)~0)'" - 2(9.36345) £ (z)' - zE(z) ps (2.3) 
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This can be seen to be related to the generalized Airy equation with a deformed 
matrix potential. 

Reversing the process of the previous subsection which was (1) matrix potential 
(2) master matrix (3) Orthogonal polynomial (4) contour integral (4) generalized 
Airy function (5) generalized Airy differential equation, one can attempt to recon- 
struct a master matrix. 

Rescaling the argument of we define: 



-(v^.gssge)) 1 / 4 ^ 

So that one has the following approximate equation for small z: 

E*(z)"' — siE*(z)' — zE*(z) ~ 

where: 

9.36345 . d . 

si = . 2.4 
V5.95896 

This appears related to the deformed (3, 1) minimal model discussed in [9] with 
deformation parameter s\, in the same way that the function ^ was related to $ in 
the undeformed (3, 1) model. 

The solution to the equation for H* is denoted by *&(z, s\) and is: 



e 4" 



We plot this function in Figure 8 and note the qualitative similarity to the 
function for small z. 
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Figure 8: Plot of the logarithm of the magnitude of the *$>(x/a, s%) function on the 
real axis for si = (5^95395)) 1/2 anc ^ a = \/2(5. 95896) 1//4 . The zeros are symmetrically 
located on the positive and negative real axis. The position of the zeros are given 
by the location of spikes pointing down. 
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One can improve the approximate equation (2.3) by including higher order terms 
in the <p expansion of U ((/)). Keeping terms up to 6 in (2.2) one obtains the 
approximate differential equation : 

6(2.15104)3"'" (z) + 4(5.95896)- (z)"' - 2(9.36345)H(z)' - zE(z) « 
Now rescaling can put the equation in the form: 

S**(z) + S^^^^^z) — Si^**(z) — z^,** (z) ~ 

with deformation parameters s\ and S3. This equation appears related to the de- 
formation of the (5, 1) minimal model of [9]. 

Finally we can define a function $>(z, si) as the solution to: 

$(z, Si y" - 8l $(z, Sl )' + z$(z, ai) = (2.5) 

which is real on the real axis and decays non-oscillatory for large positive y. Using 
the results of [9] <E>(z, s±) is the FZZT partition function associated with the matrix 
potential: 

= i^( 3 (- 1+ ^) + |(- 1+ ^) 2+ ^- 1+ ^) 3+ ^(- 1+ ^)) 

After rescaling and shifting the point of origin of the potential one can define poly- 
nomials for the matrix model deformed by the parameter si through: 

(1 -\- g i \ ^ 
— j^) 5f(exp(- r ^(3x + \x 2 + \x* + ^=x - xz))\ x=0 

The master matrix which has this as characteristic polynomial is a simple rescaling 
of the coupling constant of the master matrix of the (3, 1) minimal model and is 
given by: 

M hj = (i - l)(i - 2)5 hj+2 + 3(i - l)6 hj+1 + 1(1 + ^=)<*i+U 

This master matrix can develop complex eigenvalues for large enough N and si. In 
particular for N > 34 and s\ given by (2.4) the eigenvalues are complex. However 
the function ^f(z, s±) obtained from changing the sign of z in the third term in (2.5) 
is very different from <&(z, si) in this respect. It would be of interest to determine 
the master matrix associated with ^(z,s\) and it's corrections for terms involving 
S3 and higher, which should in principle converge to the Ricmann H function. 
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3 Other related approaches 



In this section we compare our approach with other approaches to the Riemann E 
function. 

3.1 Expansion in Meixner-Pollaczek polynomials 

In [23] the S function was expanded in Meixner-Pollaczek polynomials. These can 
be expressed as : 



. ,(2n + l)!!. n 3 .1 3 



(3.1) 



These polynomials are the characteristic polynomial of a matrix with nonzero com- 
ponents: 

Mij = i(i + -)5i, j+ i + Si+xj 
For N = 8 this matrix is given by: 
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The characteristic polynomial of this matrix is: 



363825 74247 



16 



10493 



z A - 154z 6 + / 



which agrees with (3.1) for N = 8. 

The expansion of the S function with an exponential factor can be expanded in 
terms of the Meixner-Pollaczek polynomials as [23]: 

oo 

E(z)e-™/ 4 = J2 b nPn{z) 
n=0 



Terminating this series at N one can write this expansion as the characteristic 
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polynomial of a N x N matrix. For N = 8 this is given by: 



/ 1 \ 

3/2 1 

5 1 

21/2 1 

18 1 

55/2 1 

36 1 



V b /b 8 h/h b 2 /b 8 b 3 /b 8 b 4 /b 8 b 5 /b s 105/2 + b 6 /b 8 b 7 /b 8 J 

When b n are taken to zero this reproduces the matrix (3.2). 
The coefficients b n are linearly related to the integrals [23]: 

Unlike the a n of the previous section there are closed form expressions for these 
integrals [23] . There is also some indication that there is some numerical advantage 
to the computation of b n at large n using the asymptotic expansion of the analytic 
expression. If one thinks of the origin of the S function as coming from a quantum 
mechanical system then the expansion of the function in terms of different polyno- 
mials is similar to different choices of basis functions for the quantum description. 
Although physically there is no difference for quantum physics from the choice of 
basis, numerically there is some advantage if overlap integrals can be performed 
analytically. 



3.2 Riemann-Hilbert formulation 

Besides the large N approach to matrix models one can develop a Riemann-Hilbert 
formulation of these theories [24] [25] [26] [27] [28] [29]. In [24] it is shown that one 
can use contour integrals to produce Riemann-Hilbert representations of special 
functions like the Airy function and Riemann zeta function. The jump matrix is an 
upper triangular matrix of the form: 

(J 2 "f z) ) 

and the function of interest is given by the contour integral: 

J ( z ) = / 9(<P,z)dtp 
Jc 
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or alternatively as the unique solution to the Riemann-Hilbert problem determined 
by the pair (G, C). For the Airy function g(<fi, z) is given by 

g(<p, z) = e ^ 3 / 3 +^ 

while for the Riemann zeta function it is given by 

g(<p,z) = e - u ^ +i ^ z 

Thus in the Riemann-Hilbert formulation the relation between the jump matrix 
for the Airy function and the zeta function is the same as in the previous section, 
namely the replacement of i^- with —U (</>). 

4 Conclusion 

We have discussed the master matrix formulation of FZZT partition functions. We 
derived the master matrix associated with the (2, 1) minimal model and found agree- 
ment with the results of [16]. The characteristic polynomial of the master matrix 
was the FZZT partition function of the (2, 1) model which is the Airy function 
after taking the large N limit. We also derived the master matrix of he (3, 1) mini- 
mal model and related it's characteristic polynomial with the (3, 1) FZZT partition 
function which is a generalized Airy function $>(z). In both cases the zeros of the 
FZZT partition function were on the real axis. The extension to the general (p, 1) 
minimal model should be straightforward. Interpreting the Riemann E function as 
a FZZT partition function we developed a controlled expansion of the function in 
terms of parameters si, S3, . . . of a Konsevich type potential U{4>). We showed how 
this is related to the FZZT partition function of a minimal model with deformation 
parameters given by Sj. For order s\ one uses a deformed (3, 1) matrix model to 
describe the theory. To order S3 one uses a deformed (5, 1) model etc. This proce- 
dure should converge to the Riemann S function as one includes higher and higher 
order terms in the expansion of U(<j>). More work relating the master matrix of the 
FZZT partition functions &(z, si, S3, . . .) to the function ^(z, si, S3, . . .) is needed 
to obtain quantitative and qualitative insight into the arrangement of the zeros of 
the W(z, si, S3, . . .) and H(z)in the complex plane. However a comparison between 
Figure 6 and Figure 8 is encouraging. Finally we compared our approach to other 
approaches to the Riemann E function which involve large N matrices associated 
with expansion in Meixner-Pollaczek polynomials and the definition of the Airy and 
E functions as the solution of a Riemann-Hilbert problem. 
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